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Abstract 

A supersymmetric standard model with two Higgs doublets and two Higgs singlets is 
investigated if it can accommodate the possibility of spontaneous CP violation. Assuming 
the degeneracy of the scalar quark masses of the third generation, we find that spontaneous 
CP violation in the Higgs secor is viable in our model. In the case of spontaneous CP 
violation, the masses of the lightest two neutral Higgs bosons are estimated to be 80 and 125 
GeV for some parameter values in our model, which are consistent with LEP2 data. 



I. INTRODUCTION 



The violation of CP symmetry can be explained in the standard model [1] in terms of a complex 
phase in the Cabibbo-Kobayashi-Maskawa (CKM) matrix [2] for the mixing of the quark weak 

cigcnstatcs. Meanwhile, in principle CP symmetry may be violated through the mixing between 
the scalar Higgs boson and the pseudoscalar one in the neutral Higgs sector of the models 
with at least two Higgs doublets [3] . In practice in the minimal supersymmetric extension of the 
standard model (MSSM) [4], which has just two Higgs doublets, neither explicit nor spontaneous 
CP violation is found to be possible at the tree level, because the complex phase can always be 
eliminated by rotating the Higgs fields. 

In the MSSM, explicit CP violation is viable through the radiatively corrected Higgs potential 
[5] . In explicit CP violation scenario, the presence of CP phase at the 1-loop level in the MSSM 
leads to significant modifications to both the Higgs boson mass and the mixing among them, as 
radiative corrections due to the third generation quark and scalar quark loops become important. 
Recently, the MSSM searches for explicit CP violation were extended to radiative corrections 
due to the W boson, the charged Higgs boson, and the chargino contributions [6]. The radiative 
corrections from those loops have been estimated to be significant to the mixing among the 
neutral Higgs bosons, especially if the ratio of the vacuum expectation values (VEVs) between 
two Higgs doublets is large [6]. Also CP symmetry in the MSSM can be spontaneously broken in 
the framework of the radiatively corrected Higgs sector of the MSSM. However, the spontaneous 
CP violation is disfavored because it leads to a very light neutral Higgs boson, which has already 
been ruled out by the Higgs search at the CERN e+e" collider LEP2 [7]. 

The Higgs sector of the MSSM can be extended by adding neutral Higgs singlet fields in 
addition to two Higgs doublets [8]. The simplest extention of the MSSM is next to the MSSM 
(NMSSM). It might as well be called the (M+1)SSM [8]. The primary motivation for the 
(M+1)SSM is that this model can avoid the so-called ^u-problem in the MSSM. Thus, while 
the MSSM introduces by hand the dimensional /j, parameter, the (M+1)SSM can generate the 
corresponding quantity dynamically by means of the VEV of the neutral Higgs singlet field. 

For CP violation, it has been found that, unlike the MSSM, the (M-l-l)SSM may at the 
tree level violate CP symmetry explicitly. By assuming the degeneracy of the scalar top quark 
masses in the Higgs sector of the (M+1)SSM, it is found that large explicit CP violation may be 
realized as the VEV of the neutral Higgs singlet approaches to the electroweak scale [9] . At the 
1-loop level, in the (M-l-l)SSM, the effects of explicit CP violation to the neutral and charged 
Higgs boson masses are investigated by taking radiative corrections due to the quarks and scalar 
quarks of the third generation into account [10]. We have elsewhere extended by including the 
contributions of the loops of the W boson, the charged Higgs boson, and the charginos to the 
neutral Higgs boson masses as well as the mixing among them, where the production cross 
sections for the neutral Higgs bosons are calculated via the Higgs-strahlung process at future 
e'^e~ linear collisions when the lightest neutral Higgs boson is dominantly composed of the 
singlet field [11]. 

Meanwhile, it has been discovered that spontaneous CP violation cannot be invoked at the 
tree level in the (M+1)SSM, if the stability of the vacuum of the (M+1)SSM should be ensured 
[12]. This situation has not been improved by including radiative corrections arising from the 
degenerate scalar top quark masses [13]. More precisely, in this case, the CP violating minima 
are possible only when A is very small for tan (5=1, and thus two neutral Higgs bosons becomes 
very light. In the Higgs sector of the (M+1)SSM, spontaneous CP violation may be realized by 
virtue of radiative corrections when the mass splitting effects between the scalar top quarks are 
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taken into account [14]. 

Summing up, it has been discovered that spontaneous CP violation is not ahowcd at the 
tree level for either the MSSM or the (M+1)SSM. Furthermore, these situations have not been 
improved at the 1-loop level by taking into account the contributions of the degenerate scalar 
top quark masses. Quite naturally, we would like to extend the models further to see if CP 
symmetry might be violated spontaneously at the tree level or at the 1-loop level. Hence we add 
one more neutral Higgs singlet field to the (M-l-l)SSM, to build a model of two Higgs doublets 
and two Higgs singlets. It may be called the (M-|-2)SSM. 

In this paper, we investigate whether CP symmetry can be spontaneously broken in the 
Higgs potential of the (M+2)SSM. To explore the implication of spontaneous CP violation on 
the (M+2)SSM in more detail, we concentrate on the case in which the scalar quarks of the 
third generation become degenerate in their masses. We investigate whether spontaneous CP 
violation can be realized for the Higgs sector of the (M-|-2)SSM in the case with radiative effects 
of the degenerate scalar quark masses. We would not consider explicit CP violation scenario in 
the Higgs potential of the (M-|-2)SSM. Thus, we assume that all the relevant parameters in its 
Higgs potential are real. On the other hand, in order to realize spontaneous CP violation in the 
Higgs potential, the VEVs of the two Higgs doublets as well as of the two Higgs singlets are 
assumed to be complex. 

We describe the tree- level Higgs potential of the (M-|-2)SSM in spontaneous CP violation 
scenario. Since in supersymmetric models the contribution of radiative corrections is important 
for studying the Higgs sector, we calculate the neutral Higgs boson masses as well as the mixing 
among them by using the 1-loop effective potential with the contributions of the degenerate 
scalar quark masses in the context of the (M-|-2)SSM. We find that spontaneous CP violation 
is possible for the wide parameter space of the considered Higgs potential of the (M-|-2)SSM 
without contradicting the phenomenological constraints set by the experimental data of LEP2. 

Our paper is organized as follows. In the next section, we present the Higgs potential of 
the (M-|-2)SSM at the tree level and the 1-loop effective potential with the contribution of the 
degenerate scalar quark masses. In the third section, we calculate the (M-|-2)SSM Higgs boson 
masses in spontaneous CP violation sceanrio by considering the negative experimental results 
on the Higgs search at LEP2. Conclusions are given in the last section. Appexdix supplies the 
formulae for the elements of the neutral Higgs boson mass matrix. 



II. THE (M+2)SSM HIGGS POTENTIAL 

The Higgs sector of the (M+2)SSM contains two Higgs doublets and two Higgs singlets. Let us 
introduce them: Hf = {H^, H^) and if J = (-^2^, -^2) ^■re the Higgs doublets, and A^i and N2 
are two Higgs singlets. They will yield seven neutral Higgs bosons. The most general form of the 
Higgs superpotential of the (M-|-2)SSM, containing the term with the quark Yukawa couplings 
for the third generation, is given by 

W = hbQHib% + htQH2t% + XiHleH2Ni + X2HleH2N2 

3 1 32 31^ 

where e is the usual antisymmetric 2x2 matrix with €12 = —£21 = 1) and Aj and kij = 1, 
2) are dimensionless coupling constants. The superfield = {tLjbi) consists of the left-handed 
quarks of the third generation. The superfields and 6^ denote the charge conjugate of the 
right-handed quarks of the third generation. Note that all quark and lepton Yukawa couplings 
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except for those of the third generation quarks are neglected. The presence of the cubic terms 
of the Higgs singlets in the above superpotential breaks a global U(l) Peccei-Quinn symmetry 
which leads to a pseudo-Goldstone boson at the tree level [15]. 

Prom the above superpotential, the relevant tree-level Higgs potential in the (M+2)SSM is 
obtained as 

= Vd + Vf + Vs , (1) 
where D terms, F terms, and soft terms are given as follows: 

Vf = \XiNi + X2N2fi\Hi\'^ + \H2f) 

+ \X1HfeH2 - fciiTVf - ^k22Ni - ^k2lN-,N2f 

+ \X2HfeH2 - ^kiiNf - k22Nl - ^fci27Vl7V2p , (2) 

Vs = m]j^\Hi\^ + m]i^\H2\^ + ml^\Ni\^ + ml^\N2\^ 
- {XiAx.HfeHiNi + AsA;,,/^^ e/^a^Vs + H.c.) 

'fell - ^22 , k^ 

Ai 
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+ ^Afc,,7V| + ^Ak,,mNl + ^Ak,,NfN2 + H.c.) 





where gi and 52 are the U(l) and SU(2) gauge coupling costants, respectively, a = {a^,a'^,a^) 
are the Pauli matrices, mj^. and m%. (i = 1, 2) in the soft terms are the soft supersymmetry 
(SUSY) breaking masses, and ^a, and Ak^. = 1, 2) are the trilinear soft SUSY breaking 
parameters with mass dimension. In supersymmetric models, the most important contributions 
to the 1-loop effective potential come from the loops of the quarks and scalar quarks of the third 
generation. Neglecting the mixing in the scalar quark mass matrix [13, 16], the 1-loop effective 
potential is given as [17] 

(ftglflil' + mlf log _ ftllfl.l^log I 1:11^1- I } , (3) 

where A is the renormalization scale in the modified minimal substraction (MS) scheme and the 

soft SUSY breaking masses satisfy mf = m? = m? » m? and m? = m? = m? S> m?. Thus 

tL tR t bL bR b 

the full Higgs potential V is composed of V^' and V^. 

In this model, the charged Higgs sector automatically conserves CP. Under reasonable as- 
sumptions, either explicit CP violation or spontaneous one can occur in the neutral Higgs sector 
through the mixing of the scalar and pseudoscalar Higgs bosons. We prohibit the possibility 
of explicit CP violation in the above Higgs sector. Thus all the dimensionless coupling con- 
stants and the trilinear soft SUSY breaking parameters in the Higgs potential are assumed to be 
real. For simplicity, we assume here that some of the parameters satisfy the following relations: 
Ai = A2, Ax^ = Ax^, ku = k22, ki2 = k2i, A^^^ = A^^^ and Ak^^ = Ak^^. Thus at the tree 
level rn'j^_^ and rnf^^ have same values. Of course the two soft SUSY breaking masses can take as 
different values. At the 1-loop level, the values of mj^^ deviate from those of rnj^^ by considering 
radiative effects of the degenerate scalar quark masses. 

Let us first consider the case of CP conservation. In this case, the Higgs potential should be 
such that the VEVs of the four neutral Higgs fields are ah real: {0\H^\0) = vi, (0|i?^|0) = ^2, 
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(0|A''i|0) = xi, and (0|A^2|0) = X2, after the electroweak symmetry breaking. As is well known, 
the charged Higgs fields do not develop the VEVs because the U(l) gauge symmetry has to 
be guaranteed in the vacuum. The four neutral scalar Higgs bosons can be distinguished from 
the three neutral pseudoscalar ones, as there is no mixing between them in the case of CP 
conservation. Among the four neutral scalar Higgs bosons, the one that has the smallest mass 
is referred to as the lightest scalar Higgs boson. The upper bound on its mass is expressed as 

m| cos 2/?2 + v'^{Xl + Ai) sin^ 2/3 
+ ^< log (l + d] + / 

where the ratio of the VEVs of the two Higgs doublets is given by tan (3 and the neutral gauge 



boson mass is defined as m\ = {gf + (72)^2/2 with v = \Jv\ + V2 = 175 GeV. Note that the 
above expression does not depend both on the VEVs of the neutral Higgs singlets and on the 
renormalization scale. In the above formula, the last two terms come from radiative corrections 
due to the degenerate scalar quark masses. If CP symmetry is conserved in the Higgs sector, the 
tree-level upper bound on the lightest scalar Higgs boson mass depends critically on the sizes 
of Aj = 1, 2). Notice that the above expression is essentially equal to the upper bound on 
the lightest scalar Higgs boson mass of the (M-l-l)SSM. Thus, even though an additional Higgs 
singlet is present in Higgs sector of the (M-|-2)SSM, the upper bound on the lightest scalar Higgs 
boson mass is essentially the same as the corresponding one in the (M+1)SSM. 

Now, consider the case of spontaneous CP violation in the (M+2)SSM. In order to accom- 
modate spontaneous CP violation in the Higgs sector, we have to allow complex values for the 
VEVs of the four neutral Higgs fields as 

(0|i7°|0) = vie''^' , 
(0|i/0|0) = v^e'^^ , 
(0|iVi|0) = , 
(0|iV2|0) =X2e^^* , 

where we take Vi and Xi [i = 1,2) to be real and positive. The conditions that the Higgs 
potential becomes minimum at vi, V2, xi, and X2 yield four constraints, which can eliminate the 
soft SUSY breaking masses mjj^, rnjj^, and from the potential V. 

By redefining the overall phase of the Higgs potential, a combination of the four complex 
phases can be absorbed away, leaving three independent combinations of the phases. They may 
be chosen to be 

Note that there is no relative phase among the VEVs in V^. As the Higgs potential at the 
tree level evokes the spontaneous CP violation, the vacuum of the model has to be defined as 

a stationary point with respect to the three physical CP violating phases, ^, ^1, and ^2- The 
minimum condition for the vacuum with respect to (j) yields 

SAp^jXi sin - -y-^ + d>ki2X\X2 sin(2^ — 4>i — 4>2) + ^A\iX2 sin - 

fell = X ^ ^ (j5) 

4^2 sin \ 6 — (f)i ^ j -|- yt^ ^ ^ 
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while the remaining two minimum conditions for a vacuum with respect to and <^2 niay be 
arranged to yield expressions for A^^^ and A^^^ as 

2A;iia:J sm(0i) V 3 / 3xjsm((^i) V 3 



and 



where 



2A;i2f^Ai.X2 sin2/3 . , , ^ 8/^12X2 . f (pi 



^ ^ 2 '■ / , \ — sin(2(?!) - 01 - 02) + _ 2 ■ >^ . sin 



'->2 



/ciiXi sin(0i) 9xisin(0i) V 3 

Aki2h2X2 . /01 + 202\ , V'^XIX2 . /01-02 



sm + o— — -— — sm 



3A:iia;^ sin(0i) V 3 / /ciixfsin(0i 

4kuxl . /20i-202\ , 2i;2^isin2/3 . /, , h 
■ sm H . sm — 0i — 



3xisin(0i) V 3 / xisin(0i) V 3 

2ki2Aki^X2 . /20i + 02\ , 8A;i2X2 . /0i-02\ 

■ sm + sm , (6) 



3A;iia;i sin(0i) V 3 / 9sin(0i) 



An 

A^r. = ^ , (7) 



.d _ ki2xl / 0i + 202 \ ki2xl (2(t)i + . 



3A;iia;i sin(0i) \ 3 / 8^11X2 sin(02) V 3 

2A:i2X2 . /20i + 02\ 2A;i2Xi 

— — - sm ; — - su 

3A;iiXi sin(0i) V 3 / 3kiiX2 sin(02) 



.„ 'u''AiX2 sin2/3 f . / 0i \ 

= ^7 3 • f/^ 1 4fciiX2sm - - 02 + 3^. 

6A;iixf sm(0i) I V 3 7 

f^AiXisin2/3 [ • ; (p2 



ai 3 ■ /J, \ ) 4A;iiXi sin - 01 - — + 3^Ai sin - — 
6fciiX2 sm(02) t \ 3 / V 3 

~ ni ^ 2^^^^-! ^ (s^^^Ai sin 2/3 sin(20 - 0i - 02) - 4fciiXi sin f ^^-^-^ 
9fciiX2sm(02) L V 3 

+ ^ ( 2/3 sin(20 - 0i - 02) + 4fciix^ sin 

9fciiXiSm(0i) L V 3 

~ o ■ / , X 1 Ai sm 2/3 sm - — - 02 - 2feiiXi sm 

3x2sm(02) I V 3 / V 3 

+ , , |3.^Ai sin 2^ sin f - 0i - 2fciixi sin ^^'^^ " ^'^^ 

3x1 sm(0i) L V 3 / 

V^\IX2 . /01-02\ , f^Afxi . /01-02 

+ o . . , N sm + o . . , N sm ' 



/ciixl sini 


(<^2) 


X2 1 




sin(0i) J 





/ciixf sin(0i) V 3 

+ -ki2 sm <^ . . + 

9 V 3 y lsm(02) 

In order words, the three minimum equations with respect to 0, 0i, and 02 determine ku, A/.-^^^ 
and ^fei2- They are thus no longer independent parameters but dependent on others as given 
in the above expressions, determined by the solutions of the above three minimum equations. 
Consequently, the Higgs sector of the (M+2)SSM has as free parameters 0, 0i, 02, tan/3, xi, 
X2, ^Ai, Ai, A;i2, m?, and m?. 

In order to calculate the Higgs boson masses, let us redefine the Higgs fields of the (M+2)SSM 
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in unitary gauge as 
Hi = 



(vi + Si +isuil3A) 
sin^C+* 



Ho = 



cos I3C+ 

gi(</.-,^l/3-02/3) + S2 + i COS 13 A) 



(8) 



ATi = [ e^'^i/3(xi+Xi+iyi) ] , 



No 



X2 + X2 + iY2] 



where is the charged Higgs field, and in the absence of spontaneous CP violation, Si and 
Xj (i = 1, 2) are the neutral scalar fields and A and (i = 1, 2) are the neutral pseudoscalar 
fields. The masses of the neutral Higgs bosons are obtained as the eigenvalues of the 7x7 mass 
matrix in the basis of {Si, S2, A, Xi, X2, Yi, Y2), whereas the mass of the charged Higgs boson 
is rather straightforwardly obtained by differentiating twice the Higgs potential with respect to 
the charged Higgs field. Since the Higgs sector of the (M+2)SSM is larger than that of the 
(M+1)SSM by a neutral Higgs singlet, the charged Higgs sector of the (M+2)SSM is equal to 
that of the (M+1)SSM. 

The mass of the charged Higgs boson is thus obtained as 



2 \2 2 , 



2\\A 



8Xik 



sin 2/? 



xi cos I (j) 



Y ] + X2 COS 
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+ 



3 sin 2/3 
8AiA;i2XiX2 
3 sin 2/3 



cos ^(j 



-91- 



+ X2 COS ( (p — 



cos(20 -(pi- 4>2) , 



(9) 



where the charged gauge boson mass is defined as m|j^ = 52^^/2. The tree-level mass of the 
charged Higgs boson is not changed by taking into account radiative corrections due to the 
degeneracy of scalar quark masses. Note that complex phases enter into the charged Higgs boson 
mass. This is remarkable because spontaneous CP violation does not occur in the charged Higgs 
sector but in the neutral Higgs sector. As mentioned before, the charged Higgs sector has an 
intact CP symmetry since the pair of the charged Higgs bosons are CP-even. Hence, we might 
say that the effect of the spontaneous CP violation in the neutral Higgs sector is transferred or 
shifted to the charged Higgs sector. The shift effect on the mass of the charged Higgs boson 
would vanish if CP symmetry is conserved in the neutral Higgs sector, i.e., ii (p = (pi = (p2 = 0. 
Then we have 



m 



H+ 



m 



w 



Ay + 



2XiAx-^{xi + X2) 8Xikii{xl + X2) 8AiA;i2a;iX2 



sin 2/3 



+ 



3 sin 2/3 



+ 



3 sin 2/3 



which may further be reduced to the tree-level mass of the charged Higgs boson in the (M-l-l)SSM 
with CP conservation, if the relevant parameters are set as A = Ai, = k = kn = ki2, 
and x/2 = xi = X2. 

Now, let us derive expressions for the 7x7 symmetric neutral Higgs mass matrix in the (^i, 
S2, A, Xi, X2, Yi, l2)-basis, It is calculated by differentiating twice the Higgs potential with 
respect to the seven neutral Higgs fields. The elements of the neutral Higgs boson mass matrix, 
Mij (i, j = 1 to 7) at the tree level are given in Appendix. By considering the radiative effect of 
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the degenerate scalar quark masses, the only modification to the tree-level mass matrix for the 
neutral Higgs boson is 



In the CP conserving limit, when (p = (pi = (p2 = 0, the mixing elements among the scalar 
components and the pseudoscalar ones in the neutral Higgs boson mass matrix, Mjj (i = 1, 2, 4, 
5, and j = 3, 6, 7) are zero. Then, Mij = 1, 2, 4, 5) become the elements of the neutral scalar 
Higgs boson mass matrix while = 3, 6, 7) become those of the neutral pseudoscalar Higgs 

boson one. In general, because the complex phases are nonzero in spontaneous CP violation 
scenario, Mij (i = 1, 2, 4, 5, and j = 3, 6, 7) are nonvanishing. 

Consider, for example, Afia and M23. M13 describes the mixing between and A. Similarly, 
M23 describes the mixing between S2 and A. We can reexpress M13 by substituting the minimum 
condition for kn into it as 
g 

— -^Mki2XiX2 sin{2(f> - <pi - ^2) cos/3 , 

and similarly for M23 by exchanging cos/3 ^ sin/3. One can easily notice that the source of 
nonzero M13 or M23 is —AXiHJ eH2NlN2 + H.c. terms in Vp of the tree-level Higgs potential. 
Therefore, the situation is different from the (M-l-l)SSM, where the scalar-pseudoscalar mixings 
do not occur in any of the two Higgs doublets at the tree-level (M-l-l)SSM. We would like to 
emphasize that in the (M-|-2)SSM there exists scalar-pseudoscalar mixings even at the tree level 
by the minimum conditions. In other words, if the Higgs potential of the (M-|-2)SSM at the tree 
level is such that the VEVs of four neutral Higgs fields are complex, the mixings between scalar 
and pseudoscalar Higgs fields would take place generally, which implies the occurrence of the 
spontaneous CP violation. 



Now, we have to check if our model is phenomenologically acceptable. That is, according to either 
or not the model possess a significantly large region in its parameter space that is consistent 
with the present phenomenology, the model becomes either physical or merely mathematical, 
respectively. For this purpose, we would like to examine if our model predicts a reasonable mass 
for the lightest neutral Higgs boson, and if it passes the experimental criterion set by LEP2 
data. These jobs are performed numerically by using Monte Carlo method. 

The free parameters we have in the tree-level Higgs sector of the (M-|-2)SSM with spontaneous 
CP violation are (j), 0i, 02, tan/3, xi, X2, '^Ij ^12, and m?. For our numerical analysis, 
the ranges for the parameters are set as follows: < (/>, (/>i , (/>2 < tt, 2 < tan/3 < 40, < 
xi, X2, A\-^ < 500 GeV, < Ai < 0.87, and < ki2 < 0.63. These ranges define the extent of 
the parameter space we will explore. We fix m| and m? in the effective potential at 1000 GeV. 

Let us first estimate the neutral Higgs boson masses by diagonalizing the 7x7 symmetric 
mass matrix of the neutral Higgs bosons. The seven eigenvalues are sorted such that nih- < rrih. 
for i < j. We scan 10^ points in the relevant paramter space and numerically evaluate nih. [i 
=1 to 7). For each point, we also evaluate the mass of charged Higgs boson. We find that the 




(10) 



III. NUMERICAL ANALYSIS 
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masses of neutral Higgs bosons as well as the charged Higgs boson mass are positive for the 
most part of the parameter space of the (M+2)SSM that we consider. Consequently, we note 
that no serious reduction occurs in the considered parameter space by either the spontaneous 
electroweak symmetry breaking or the minimum conditions for the vacuum. 

Next, we consider the production of a neutral Higgs boson in e'^e~ collisions. For the 
center of mass energy of LEP2, the neutral Higgs bosons are dominantly produced through the 
Higgs-strahlung process e+e~ — > — > Zhi (i = 1 to 7) and the Higgs-pair production process 
e^e" Z* ^ hihj {i, j = 1 to 7, where i j)- The relevant coupling coefficients for the these 
processes are 

2 2 

Glzh, = j^^^(cos/30a + sin/30i2)' , 
' 4 cos^ Uw 

G\h,h, = , %n [{Oj3(cos I30i2 - sinpOn) - Ois{cos (30j2 - sm/30ji)}^ 

■' 4 COS uw 

+ {OnOji + Oi20j2 + OisOi^}^] , (11) 

where 9w is the Weinberg mixing angle, and Oij is the elements of the orthogonal transformation 
matrix that diagonalizes the mass matrix of the neutral Higgs bosons. Note that the expression 
for the coupling coefficient in the Higgs pair production is different from that of the MSSM. The 
production cross sections of the Higgs bosons receive the effects of spontaneous CP violation 
and the radiative effects of the degenerate scalar top quark masses through the above coupling 
coefficients via Oij. 

Although candidates of detecting a Higgs boson has been suggested by LEP2 [18], no concrete 
evidence is available yet; thus, the parameter space of the (M+2)SSM may be constrained 
experimentally if we assume the discovery limit of 0.1 pb at LEP 2 with y/s = 200 GeV. Let us 
denote at the total cross section for the neutral Higgs productions via the above two processes. 
The total cross section is evidently dependent on the neutral Higgs boson masses. For a given 
set of parameter values, we evaluate m/j. numerically in order to insert into the formula of at- 
Then we evaluate at numerically. We carry out the analysis for the whole parameter space we 
have defined, by choosing randomly a total of 10^ points in the space. Finally, we select those 
points which yield at < 0.1 pb. They are consistent with experimental constraints and thus 
establish a subset of parameter space that are phenomenologically acceptable. 

We find that there is an indirect lower limit of 110 GeV for the mass of charged Higgs boson 
by considering in this way the failure of the neutral Higgs boson production at LEP2 with ^/s 
= 200 GeV. This might be considered as its phenenoniclogical lower bound in our model. Since 
the two neutral Higgs singlets in our model distinguish it most distinctively from (M+1)SSM, 
it is worthwhile to study the behaviors of their VEVs. In Fig. 1, we show the points that are 
accepted by the above procedure in [xi, X2) plane. It is interesting to notice in Fig. 1 that both 
xi and X2 cannot simultaneously be small: Roughly, they cannot simultaneously be smaller than 
- 20 GeV. 

In Fig. 2(a), we plot (solid curve) and m/jj (dashed curve) as a function of tan/3, for 
^ = 7r/10, = 27r/3, (/)2 = 7r/20, xi = 200 GeV, X2 = 400 GeV, Ax, = 250 GeV, Ai = 0.1, 
ki2 = —0.5, and = = 1000 GeV. In the case of Fig. 2(a), the mass of the charged Higgs 
boson is estimated to be between 287 and 962 GeV and at between 0.1 and 0.4 fb as tan/3 goes 
from 2 to 40. In Fig 2(b), we repeat the same calculation as Fig 2(a), except for xi = 400 GeV 
and X2 = 200 GeV. In the case of Fig. 2(b), we have 230 < mH+ < 755 GeV and 8.6 < a* < 53.6 
fb for 2 < tan/3 < 40. For larger tan/3, at becomes smaller and larger. Those parameter 

values are consistent with the phenomenological constraints of LEP2 at s/s = 200 GeV. Figures 
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2 display that two neutral Higgs bosons /ij (i = 1, 2) are relatively heavy for the allowed range 
of 2 < tan/3 < 40. Those results are quite dissimilar to the case of the (M+1)SSM. In the case 
of the (M+1)SSM, the sum of two lightest neutral Higgs boson masses is smaller than 100 GeV 
and the charged Higgs boson mass is smaller than 110 GeV [13]. 

IV. CONCLUSIONS 

It is known that in super symmetric models with at least two Higgs doublets CP symmetry may 
be violated through the mixing between the scalar Higgs boson and the pseudoscalar one in 
the neutral Higgs sector. For the trcc-lcvel Higgs potential of the (M+1)SSM, spontaneous CP 
violation cannot occur in its neutral Higgs sector since the CP violating vacuum is not stable 
and thus the mass of the lightest neutral Higgs boson is negative for a reasonable part of the 
parameter space. 

Within the context of the (M+1)SSM, spontaneous CP violation scenario is hardly improved 
when the 1-loop effective potential due to the degenerate scalar top quark masses is taken into 
account. We have investigated whether or not the above situation might be effectively improved 
by considering the (M+2)SSM with the 1-loop effective potential due to the degenerate scalar 
quark masses for the third generation. The most general form of the tree-level Higgs potential 
is constructed, where the interaction terms between the two Higgs singlet fields are included. 
We also take into account the radiative effects due to the degenerate scalar quark masses of the 
third generation to the tree-level Higgs potential. 

The Higgs potential is supposed to develop complex VEVs for the four neutral components 
of the Higgs multiplets, in order to accommodate spontaneous CP violation. From the complex 
VEVs, three complex phases appear as independent free parameters of the model. With respect 
to these complex phases, three minimum conditions for the Higgs potential are calculated to 
reduce the number of free parameters by three. Then, the mass matrix for the seven neutral Higgs 
bosons are calculated. In the (M+2)SSM, the matrix elements (M13, M23) for Si- A mixing and 
S2-A one would not vanish for spontaneous CP violation scenario. The presence of the mixing 
elements is found to be essential for spontaneous CP violation. 

By assuming the discovery limit of 0.1 pb at LEP2 with ^/s = 200 GeV, we explore the 
parameter space of the Higgs sector of the (M+2)SSM, if the lightest neutral Higgs boson might 
have escaped from detection at LEP2. Our model is numerically analyzed if it is phenomenolog- 
ically acceptable. In spontaneous CP violation scenario, the mass of the lightest neutral Higgs 
boson is evaluated to be as large as about 80 GeV for some parameter space. The charged Higgs 
boson has 110 GeV as the phenomcnological lower bound on its mass in our model. Conclud- 
ingly, we have found that spontaneous CP violation in the (M-|-2)SSM can be accommodated 
for the Higgs potential including radiative effects due to the degenerate scalar quark masses. 
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APPENDIX: HIGGS BOSON MASS MATRIX 

At the tree level the elements of the neutral Higgs boson mass matrix are expressed as 
follows: 

Mil = (mzcos^)^ + Ai^Ai ja^icos ^</) - + X2 cos - I tan/? , 



+ ^Ai^ii jxi cos (^(p- 4>i- y j + xl cos (^<?^ ~ Y ~ '^2j j tan /3 
4 

+ o '^1^122^13^2 cos(20 - 01 — <p2) tan/3 , 
Mi2 = 2i;2A? sin 2^ - ^ sin 2/3 - Ai^Ai |^i cos " y) + cos (-^ " y) } 



— ^Ai/cii jx^cos 1^0- 01 - +^2COS ^ 
4 

— 2 AiA;i2Xia;2 cos(2(/) - 01 - 02) , 

Mi3 = Ai^Ai ja^i sin (^(p - + 2;2 sin ^0 - Y ^ I cos P 

4 (of 02\ 2 / 01 M 

+ -^Mku jcci sm I - 01 - — I + X2sm ( - — - 02 I j- cos/3 
4 

+ ■^■^1^122^13^2 sin(20 — 01 — 02) cos/3 , 
Mi4 = 2vXl |xi + X2 cos — ~^ I COS /3 — vAi^Ai cos ^0 — sin/3 

— ^'wAifeiicci cos ^0 - 01 - sin/3 - ^?;Aifci2a;2 cos(20 - 0i - 02) sin/3 , 
Mi5 = 2'uAi |a;2 + xi cos (^^^~^~^^ | cos /3 — ■uAi^Ai cos ^0 — sin/3 

— -f Aifciia:;2 cos f ~ ^ ~ 02^ sin /3 — -?;Aifci2a::i cos(20 — 0i — 02) sin /3 , 
Mi6 = uAiAai sin ^0 - sin/3 - ^i;AiA;iixi sin ^0 - 01 - sin/3 

— ^i'AiA;i2a;2 sin(20 — 0i — 02) sin /3 — 2vAiX2 sin cos /3 , 
Mi7 = i;Ai^Ai sin ^0 — y^ sin/3 — -wAi/ciiX2 sin ~ y" ~ "^2^ sin/3 



— -■uAiA;i2a;i sin(20 — 0i — 02) sin /3 + 2vAfxi sin ( ) cos /3 , 

M22 = Mil (cos /5 sin /3, tan /3 cot /3) , 

M23 = Mi3 (cos /3 ^ sin/3) , 

M24 = Mi4(cos/3 sin/3, sin/3 — > cos/3) , 

M25 = Mi5(cos/3 ^ sin/3, sin/3 — *• cos/3) , 

-^26 = Mi6(cos /3 ^ sin /3, sin /3 ^ cos /3) , 

M27 = Mi7(cos/3 ^ sin/3, sin /3 ^ cos /3) , 
Ai^Ai^;! 



M33 = 



4 sin 2/3 



cos (4/3 - + ^) - cos (4/3 - 20 + ^ + ^02) I 

2Ai^r / 02\ /, 01 
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3sin 2,3 H + f + + I* N } 

SAifeii 12 /^i, J. 4>l\ ,2 ( A 



SAlfell r 2 . <?!>2\ 2 

3 sm 2p 

M34 = vAiA;^^ sin ^0 - + ^uAi/ciiXi sin ^0 - 0i 

4 

+ -^vXiki2X2 sm{2<p - 01 - 02) , 

M35 = M34(0i ^ 02, 02 01, Xi ^ .T2, 2:2 ^ Xi) , 



3 '^^ 



Mse = fAi^lAi cos ^0 - yJ - -t;Aifciixi cos ^0 - 01 ^ 
4 

- -f Alfcl2.T2 COS(20 - 01 - 02) , 
M37 = M36(01 02, 02 01, Xi X2, X2 ^ Xi) , 

M44 = — A;fiXi - kuAkj^^xi cos(0i) + — fciiA;i2a;2(3,T^ - x^) cos 

^^2,, 30 70i^'02\ /0 01 02 

+ AifcnCOS^2/3-- + — + -jsm(^----- 

^%2 , t-^Ai / 02\ . 

- - A,X2 COS [^-^ j + COS ^0 - - j sm 2/? 

+ ^i|^xicos (^^1±^) + |!aiA;i2X2COs(20 - 01 - 02) sin 2/3 , 

M45 = fAiCOSI ^ I + -«;ilKl2(Xi + X2)C0S I I 

+ ^^11^1X2 COS ^^^L_?^^ _ ^■u^AiA;i2 cos(20 - 0i - 02) sin 2/3 

^16 2 2 r /20i+02\ /0i + 202\\ 
+ ^^12X1X2 - -/ci2^fci2 j a;i cos I j + X2 cos ( I | , 

M46 = 2A;ii^jtiiXi sin(0i) - ^v^Ai^n sin ^0 - 0i - -y^ sin2/3 

- ^t;,,^ sin f?*^) _ '-^kuku.,.. sin ^ * 



3 11 Z 3 J g ..... . 3 

2 . /20i + 02\ 
+ 3Ki2^fci2a;2Sin I I , 

M47 = ^2A2,in(^^L_^^+8fc^^fc^,(3^2+^2)3i^(^^ 



,2 . /201 - 202\ 2 2, 



3 

+ |A;nXiX2 sin j _ |i;''AiA;i2 sin(20 - 0i - 02) sin 2/3 

2 r . /20i + 02\ . ^01 + 202 

+ -«^i2^fci2 1 a;i sm 1 ^ ) + X2 sm 



3 V 3 J \ 3 

M55 = M44(0i 02, 02 ^ 01, Xi X2, X2 ^ Xi) , 
M56 = M47(0i ^ 02, 02 ^ 01, Xl X2, X2 ^ Xl) , 
M57 = M46(0l 02, 02 ^ 01, Xl ^ X2, X2 ^ Xl) , 

Mee = 3fcii74fe^^xi cos(0i) - - — feiifci2X2(xi + xl) cos 



9x1 ^ ^' V 3 
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8 2 2 ( 2d)i-2d)2 \ fh~h 
-fc]^;^X2COS A]^a;2Cos 



3 V 3 ; .Ti 

,2 



+ ;r — Ai^Ai COS ( — — 1 siii 2/3 + - — AiA;i2a;2 cos(2^ - 4>i - 4>2) sin 2/3 
+ ^.^Aifcn sin (2/3 -2^+1^ + ^ 

^ fci2 . r, /2^i+^\ /0i + 2^\\ 

+ ^^'Aifcn |2sin (^2/3 + <^ - - ^) + sin (^2/3 - + + ^ 

Ml - 02 



M67 = v'^Xl COS (^^-3-^) + ^^11^12(2;? + xl) COS 

8,2 f2(j,i-2(f>2\ , 2 2, 



+ -^kiiXiX2 COS ( ^ ) ^ 3^ ^^1^12 cos(2(/) - 0i - </'2) sin 2/3 



2 f /2^i+^\ /( 
+ 3 fci2-4fc,2 I Xl cos I 1 + a;2 cos I - 

M77 = M66(01 02, 02 ^ 01, Xl X2, X2 ^ Xi) . 



+ 202 
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Figure Captions 



Fig. 1 : The points in {xi, a;2)-plane that satisfy at < 0.1 pb are plotted for < X2 < 100 GeV 
by random number generation method. The ranges of other parameters are: < (f>,(f>i,4>2 < ir, 2 
< tan/? < 40, < ^Ai < 500 GeV, < Ai < 0.87, and ku < 0.63. 

Fig. 2(a) : m/j^ (sohd curve) and m/jj (dashed curve) are plotted as a function of tan/3, for 
<f> = vr/lO, (f>i = 27r/3, (f>2 = vr/20, xi = 200 GcV, X2 = 400 GeV, Ax^ = 250 GeV, Ai = 0.1, 
ki2 = —0.5, and = m-^ = 1000 GeV. The plotted results satisfy at < 0.1 pb. 

Fig. 2(b) : The same as Fig 2(a), except for xi = 400 GeV and X2 = 200 GeV. 
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Fig. 1: The points in (xi, 2;2)-plane that satisfy at < 0.1 pb are plotted for < xi, X2 < 100 GeV 

by random number generation method. The ranges of other parameters are: < 0, 0i, 02 < tt, 
2 < tan/? < 40, < Ax^ < 500 GeV, < Ai < 0.87, and ku < 0.63. 
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Fig. 2a: (solid curve) and m/jj (dashed curve) are plotted as a function of tan/?, for 

4) = vr/lO, = 27r/3, (1)2 = vr/20, xi = 200 GeV, X2 = 400 GeV, Ax^ = 250 GeV, Ai = 0.1, 
ki2 = —0.5, and = = 1000 GeV. The plotted results satisfy at < 0.1 pb. 
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Fig. 2b: The same as Fig 2(a), except for xx = 400 GeV and X2 = 200 GeV 
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